
%input is a binary image with ones where there are positive pixels
function r = getLineOrientation(im)

distanceThreshold = [25,150];   %in pixels - opposing vessel wall has to be within these distances
angThreshold = 0.17;            %in radians - opposing vessel wall has to have an absolute normal orientation within this range

orientations = deg2rad(skeletonOrientation(im,7));   %can increase the second argument to make things a bit better

